First-principles study of high-field piezoelectricity in tetragonal PbTiOs 
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We calculate from first principles the nonlinear piezoelectric response of ferroelectric PbTiOs for 
the case of a polarization-enhancing electric field applied along the tetragonal axis. We focus mainly 
on the case of fixed in-plane lattice constants, corresponding to epitaxially constrained thin films. 
We find that the dependence of the c/a ratio on electric field is almost linear in the range up to 
500 MV/m, with little saturation. This result contrasts with expectations from Landau-Devonshire 
approaches based on experimental results obtained at lower fields, but is in qualitative agreement 
with a recent experiment in which higher fields were attained using pulsed-field methods. We also 
study cases in which the in-plane epitaxy constraint is removed, or an artificial negative pressure 
is applied, or both. These calculations demonstrate that PbTiC>3 can show a strikingly non-linear 
piezoelectric response under modified elastic boundary conditions. 

PACS numbers: 77.65.Ly, 77.84.Dy, 77.65.Bn, 63.20.dk 



I. INTRODUCTION 

Piezoelectricity was discovered in the late nineteenth 
century and has been successfully described at the phe- 
nomenological macroscopic level since the early twentieth 
century. Piezoelectric materials find a variety of tech- 
nological uses including sonar detection, electromechani- 
cal actuators and resonators, and high- precision applica- 
tions such as in microbalances and scanning-probe micro- 
scopes. Most of these conventional applications depend 
on the linearity of the strain response to the applied elec- 
tric field or vice versa. Application of large strains or 
strong electric fields may, however, generate a nonlinear 
response. Such nonlinearities are less well characterized 
and understood, despite the fact that they may have im- 
portant practical implications for degrading or improving 
device performance. 

The intrinsic nonlinear piezoelectric properties of a ma- 
terial may be difficult to access experimentally for two 
reasons. First, application of large stresses or strains may 
cause cracking or other forms of mechanical failure, while 
application of strong electric fields may cause dielectric 
breakdown, leading in cither case to a limited ability 
to access the nonlinear regime. Second, most strong 
(and strongly nonlinear) piezoelectrics are ferroelectrics, 
whose spontaneous polarization is strongly modified by 
applied strains or fields. Real samples of ferroelectric 
materials usually break up into domains having different 
directions of the spontaneous polarization, and the ob- 
served piezoelectric response may be dominated by ex- 
trinsic processes associated with changes in domain pop- 
ulations under applied strain or field. Studies of single- 
domain samples arc difficult; even if such a sample is ob- 
tained, the domain may tend to reorient under applied 
strain or field. In the case of a strong applied electric 
field, one expects that the sample will eventually be poled 
into a single-domain state, so that a subsequent increase 
of the field strength will access the intrinsic piezoelectric 



response. However, one must still overcome the problems 
associated with dielectric breakdown mentioned earlier. 

In the last two decades or so, predictive quantum- 
mechanical descriptions based on first-principles 
electronic-structure calculations have been successfully 
developed and applied to piezoelectric and ferroelectric 
materials (see Ref. [l| for a recent review). These 
methods allow one to isolate the intrinsic contributions 
to the nonlinear piezoelectric behavior of a material 
by calculating the response explicitly. We apply such 
methods here to the study of tetragonal PbTiC>3 (PT) 
thin films under a strong electric field applied parallel to 
the polarization direction. 

Our work is motivated in part by recent experimental 
and modeling studies of PT and PbZr^Tii^Oa (PZT) 
thin films.— Based on experimental measurements^ on 
PT and PZT films in fields up to -50 MV/m, Landau- 
Devonshire models were developed to systematize the 
data and extrapolate to higher fields^ For pure PT, these 
papers show dropping by about 16% and 10% as 
the electric field is increased to 50 MV/m for the cases 
of free-stress and epitaxially-constrained samples respec- 
tively. (For the latter, the in-plane lattice constant was 
constrained to that of the SrTiC>3 substrate.) These re- 
sults are roughly in agreement with previous results of 
a first-principles calculation 5 for the free-stress case, al- 
though the methods used there involved approximations 
that are removed in the present work. 

More recently, Grigoriev et alA measured the piezo- 
electric response of epitaxially constrained PZT 20/80 
(i.e., x—Q.2) films using a novel approach in which di- 
electric breakdown was avoided by the use of ultrashort 
electric-field pulses in a thin-film geometry. This al- 
lowed access to much larger electric fields than previ- 
ously possible, up to about 500 MV/m compared with 
—50 MV/m studied in the earlier experiments. The strain 
response showed almost no saturation up to the highest 
fields reached in the experiment, around 500 MV/m, con- 
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trasting with the conclusions of Refs. 2,3]. These results 
may suggest that the previous work, limited as it was to 
smaller fields and longer time scales, may have been more 
sensitive to extrinsic effects such as incomplete poling of 
the domains. 

Clearly information from predictive first-principles 
calculations would be very useful here. A preliminary 
application of such methods to the case of stress-free 
PT under strong applied electric fields^ showed not only 
a lack of saturation, but even an enhancement of the 
piezoelectric response, with ^33 increasing up to fields of 
about 550MV/m. This unexpected^ enhancement was 
ascribed to a field-induced structural transition in which 
the PT crystal adopts a supertetragonal state at high 
fields, with an extreme axial c/a ratio of 1.2-1.3. Such a 
supertetragonal state was first identified theoretically by 
Tinte et al.jL who predicted its occurrence in PT under 
conditions of fictitious negative pressure. Interestingly 
a similarly strong tetragonal distortion was recently re- 
ported for other materials, e.g. BiCo03 in its bulk ground 
stated and BiFe03 films under a compressive epitaxial 
strain. 9 These results suggest that the peculiar electrome- 
chanical response of PT reported in Ref. i might be a 
rather general property of many perovskite materials. In 
particular, it is not unreasonable to think that, by under- 
standing the interplay of composition, epitaxial strain, 
pressure and applied external fields, one might be able 
to control the crossover between the competing "normal" 
and "supertetragonal" phases in a given Pb- or Bi-based 
compound. This might produce exciting and unusual ef- 
fects, such as giant piezoelectric responses and radical 
magnetoelectric couplings (for example, a magnetic mo- 
ment collapse and insulator-to-semimetal transition un- 
der pressure was recently reported for BiCoOs 10 ), with a 
wide range of implications for technology and fundamen- 
tal science. 

To help shed light on the above issues, in this work 
we extend the analysis of Ref. H to study the impact of 
the elastic boundary conditions on the electromechanical 
response of PT. In particular, we address here the case of 
epitaxially constrained or stress-free films, with or with- 
out external pressure, and in a wide range of applied 
electric fields. For the case of an epitaxially constrained 
film, our first-principles calculations indicate that ^33 of 
PT monotonically decreases with field, but only slowly: 
by about 2-3% up to 50MV/m, 11-14% at 500MV/m, 
and 40-50% at 1500 MV/m. Our results thus suggest 
that, based on intrinsic material properties, the strain 
enhancement can continue to quite large fields, with only 
a very slow saturation of ^33, at variance with previous 
expectations 2 -^ and strengthening the findings in Ref. |4|. 
In the absence of an epitaxial constraint, we confirm the 
previous results of Ref. 6], finding a maximum of ^33 
around 550 MV/m and then a gentle drop-off to about a 
third of its highest value at 1400 MV/m. We further in- 
vestigate the field-induced strain response of the system 
under the application of a fictitious negative pressure. 

In comparing our work with experimental studies, it is 



useful to keep in mind the context and implicit assump- 
tions of our calculations. Following the standard frame- 
work of density-functional theory, we treat pure defect- 
free PT in perfect tetragonal symmetry at zero tempera- 
ture, with and without in-plane epitaxial constraint. We 
preserve the primitive periodicity, so that domain forma- 
tion is excluded. We study the response of the system 
on a time scale on which lattice and strain response can 
occur, but domain dynamics cannot occur. We can go 
to much larger fields than can be accessed even by the 
pulscd-field methods, 4 limited only by the intrinsic break- 
down fields associated with our method for treating the 
electric fiekh 11 ' 12 We thus study a regime of intrinsic be- 
havior that, we believe, is closer to that of Ref. [H than 
that of Refs. HQ. 

The paper is organized as follows. We briefly describe 
the computational details in Sec.|IT] We then present the 
results of our calculations in Sec. IIHI and compare them 
with experiment. We also describe further investigations 
of the system as we go beyond the experimental condi- 
tions, and discuss the implications of the work. We then 
summarize in Sec. IIVI 



II. COMPUTATIONAL METHODS 

The calculations are performed using density- 
functional theory with two ab-initio computer code 
packages, ABINITH and LAUTREC^ We use the 
Ceperley-Alde r 15 i 16 exchange-correlation, implemented 
in the Perdew-Zunger— and Perdew-Wangi^ parameter- 
izations for ABINIT and LAUTREC respectively. For 
ABINIT we use norm-conserving pseudopotentials gen- 
erated using the method of Ramer and Rappe^ as im- 
plemented in the OPIUM package^ while in LAUTREC 
we use projector augmented-wave (PAW) potentials. 21 
In both cases, the semicore 3s and 3p orbitals of Ti, 
and the 5d orbitals of Pb are treated as valence elec- 
trons. Plane-wave cutoffs of 50 and 30 Hartree are cho- 
sen for ABINIT and LAUTREC respectively (the PAW 
potentials being softer than the norm-conserving ones). 
The Brillouin zone is sampled by a 4 x 4 x 4 Monkhorst- 
Pack 2 ^ fc-point mesh for ABINIT and a 6 x 6 x 6 mesh for 
LAUTREC. A stress threshold of 2x 10~ 2 GPa is used for 
cell relaxation^ and forces on ions are converged below 
2.5xl0" 3 eV/A. 

In ABINIT the electric polarization is calculated using 
the Berry-phase approach 2 ^ and is coupled to a fixed elec- 
tric field Ei^. In LAUTREC the electric polarization 
is computed using the centers of the "hermaphrodite" 
Wannier orbital s 26 ' 27 and the electric displacement field 
D is used as the independent electrical variable i 6 -^ In 
either case, the appropriate E or D field is applied, and 
the internal coordinates and unconstrained lattice con- 
stants are allowed to relax within the constraints of the 
enforced PJ^mm tetragonal symmetry. In both cases, the 
results are presented as a function of E; for the case of 
LAUTREC, this is determined at each D from the com- 



3 




Electric field (MV/m) 

FIG. 1: Out-of-plane strain response to applied electric field, 
relative to zero-field ferroelectric state, under in-plane epi- 
taxial constraint. Experiment: data of Grigoriev et alA for 
several PZT 20/80 capacitors of varying thickness. Theoreti- 
cal points with fitted lines: computed results for pure PbTiOs 
using two computer codes. Inset shows the ABINIT results 
over a more extended range of fields. 

puted polarization using E — D — 4nP. 

As we shall see, the results obtained using ABINIT and 
LAUTREC are generally consistent, but with some quan- 
titative differences between them. Tests indicated that 
the codes produce almost identical results when using 
the same potentials, so we can ascribe these differences 
almost entirely to the different potentials used (norm- 
conserving for ABINIT vs. PAW for LAUTREC). 

III. RESULTS 

A. Piezoelectric response with epitaxy constraint 

Both the ABINIT and LAUTREC computer codes 
were used to compute the response of the PT system 
to electric fields up to 1400 MV/m applied parallel to 
the polarization (along the tetragonal axis). The epi- 
taxy constraint was enforced by fixing the in-plane lat- 
tice constant a to remain at a value 1.14% smaller than 
the computed equilibrium ao of bulk tetragonal PT; this 
factor was chosen to facilitate comparison with the ex- 
periments of Grigoriev et alA in which the PZT 20/80 
film was compressively strained by 1.14% relative to its 
bulk equilibrium in-plane lattice constant by its epitax- 
ial coherence with the SrTiOa (ST) substrate. Having 
computed the equilibrium ao to be 3.892 and 3.854 A us- 
ing ABINIT and LAUTREC respectively, we therefore 
set the constrained a to be 3.848 and 3.811 A for the two 
codes respectively. (We have traced the shift of lattice 
constant between codes to the choice of Ti pseudopo- 
tential, but trends are well reproduced. For example, 
the respective c/a ratios for the two codes are 1.085 and 
1.083 at the constrained a, and 1.047 and 1.044 at the 



equilibrium ao.) For each value of the applied field, we 
computed the relaxed structure subject to the epitaxial 
constraints, and computed the strain 773 = (c — co)/co, 
where c is the z lattice constant at the given field and cq 
is the zero-field value. 

The main panel of Fig. Q] shows the results in the range 
up to 500 MV/m, together with a comparison with the 
experimental data of Grigoriev et alA on PZT 20/80, 
while the inset shows the computed behavior over the 
full range. The higher fields, above ^500 MV/m, are 
still outside the range that is achievable experimentally, 
even with pulsed-field techniques. We find that the 
strain increases monotonically over the entire range up to 
1400 MV/m, with only a small tendency to begin saturat- 
ing at the highest fields. In the range up to 500 MV/m, 
the results look very nearly linear. We obtain an excellent 
fit to the results with a simple quadratic form t](E) = 
aiE + (I2E 2 (where 773 = 77 and E$ = E henceforth). 
We obtain fitted values of ai ~ 0.063 - 0.071 m/GV and 
C12 — —0.01 (m/GV) 2 , where the quoted range reflects 
the choice of computer code. 

The slopes of the curves in Fig. Q] are, of course, re- 
lated to the piezoelectric response of the material. At 
finite field it is possible to define two slightly different 
piezoelectric coefficients (I33 = drj/dE and = dc/dV, 
where V = cE is the potential drop across a unit cell; 
from drj = c~ 1 dc it follows that the definitions are re- 
lated by dg 3 = dgg 1 + E. We report 033 values here, 
consistent with the conventions of Ref. 4. However, the 
two definitions coincide at E = and differ by only about 
2% at 500 MV/m, so the difference is not significant in 
what follows. The results are presented in Fig. [2ja); we 
find a monotonically decreasing trend with electric field 
for both ABINIT and LAUTREC. However, the fall-off 
is quite slow, decreasing by only about 2-3% up to a field 
of 50 MV/m. 

Other properties of PT under applied field show a sim- 
ilar, nearly linear behavior. In Fig. Hfb) we present the 
variations of the internal coordinates as a function of ap- 
plied field (here computed using ABINIT, but similar 
results are obtained with LAUTREC). As expected, Pb 
and Ti ions, being positively charged, displace along the 
direction of the applied electric field, while O atoms dis- 
place in the opposite direction, thus increasing the fer- 
roelectric mode amplitude. The Born effective charges 
Z* , shown in Fig. 02a), are almost independent of field, 
with only a mild reduction in their magnitudes with in- 
creasing field. However, at zero field the Z* values under 
the epitaxial constraint are noticeably smaller in mag- 
nitude than their free-stress counterparts. 28 Since the 
Z*'s are almost constant, we expect the electric polar- 
ization to increase with field in proportion to the dis- 
placements shown in Fig.^b), and Fig.^b) shows that 
this is indeed the case. A satisfactory fit of the form 
P{E) = c +c 1 E+c 2 E 2 is obtained with c = 0.94 C/m 2 , 
ci = 0.39 C/(GV)m and c 2 = -0.072 C/(GV) 2 . None of 
these quantities show any anomalous change with electric 
field. 
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FIG. 2: Variation of properties of tetragonal PbTiO;? with 
electric field under in-plane epitaxial constraint, (a) Piezo- 
electric coefficient (calculated with both codes), (b) Out- 
of-plane displacements of ions from zero-field positions (cal- 
culated with ABINIT). 



We now return to a discussion of the strain response. 
Our theory is in qualitative agreement with the recent 
pulse-field data, 4 which show a nearly linear behavior 
similar to that predicted theoretically. The theory has a 
somewhat larger slope; our zero-field ^33 of ~68pm/V 
can be compared with their 45pm/V. The agreement 
seems reasonable given that the materials are different 
(PT vs. PZT 20/80). Hints of a hump-like nonlinearity 
in the experimental data of Ref. around 200 MV/m are 
not supported by the theory; if such a feature is really 
present, we would argue that it must arise from extrinsic 
effects not considered by the theory. 

Our zero-field ^33 is also significantly larger than the 
value of 45pm/V obtained by the Landau-Devonshire 
theory for PZT 20/80^ which is in reasonable agreement 
with the zero-field value of Grigoriev et al. (For detailed 
comparisons, it should be kept in mind that the misfit 
of —0.5% reported in Ref. d is only about half of that in 
Ref. 0.) However, our theoretical curves and the experi- 
mental data of Grigoriev et al. both show a much slower 
saturation of the strain response with field than was 
found by the Landau-Devonshire approach, where ^33 
was predicted to fall by about 15% already at 50 MV/m 
for PZT 20/80 with epitaxy constraint. Our c?33 falls by 
only about 2-3% over the same range. The lack of satura- 
tion means that we predict large strains of about 3.0% at 
fields of around 500 MV/m, and we also predict a large 
polarization of about 1.10 C/m 2 at such fields. Possi- 
ble reasons for the discrepancies between our results and 
those of the Landau-Devonshire theory will be discussed 
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FIG. 3: Variation of properties of tetragonal PbTiOa with 
electric field under in-plane epitaxial constraint, calculated 
with ABINIT. (a) Born effective charges, (b) Polarization. 



further in Sec. [IlTCl 

The Landau-Devonshire theory^ also provided a zero- 
field value of 31 pm/V, falling by about 10% at 50 MV/m, 
for pure PT on an ST substrate. Recall that our calcu- 
lations above were carried out at an in-plane lattice con- 
stant that was reduced by 1.14% relative to that of PT in 
order to model the strain state of PZT 20/80 in coherent 
epitaxy on a ST substrate. Therefore, in order to obtain 
a more direct comparison for the pure PT case, where 
the epitaxy with the ST substrate is almost perfect, we 
repeated our calculations with the in-plane lattice con- 
stant fixed to the theoretical equilibrium lattice constant 
of PT, finding a zero-field ^33 of 50 pm/V. This is reduced 
somewhat from our value of 68 pm/ at —1.14% misfit, but 
still quite a bit larger than their value of 31 pm/V. Again, 
we find a very slow saturation (even slower that for the 
smaller in-plane lattice constant), in disagreement with 
the Landau-Devonshire theory. 



B. Application of negative pressure and removal of 
epitaxy constraint 

In this section we present the results of further investi- 
gations of the behavior of the system under modified elas- 
tic boundary conditions that are not directly relevant to 
the thin-film experiments. These investigations are mo- 
tivated in part by previous ab initio calculations of Tinte 
et air- who showed that the c/a ratio and polarization of 
PT undergo an anomalous and strongly nonlinear varia- 
tion as a function of an artificial applied negative pres- 
sure. The c/a ratio was found to increase gradually until 
the negative pressure reached about —4.8 GPa, where the 
c/a ratio rapidly increases to ^1.20, with a subsequent 
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FIG. 4: Field-induced strain, relative to zero-field zero- 
pressure state, for tetragonal PbTiOa under in-plane epitaxial 
constraint, calculated with LAUTREC. 



slower increase up to ~1.25 at — 7GPa. Neither an elec- 
tric field nor an epitaxy constraint were applied in that 
calculation, but the results nevertheless hint at a possi- 
ble supertetragonal state of PT that might be accessed 
under unusual boundary conditions. We also note that 
negative pressure is sometimes used to simulate the effect 
of chemical substitution, and that such an approach may 
be relevant here for providing hints about the behavior 
of PZT solid solutions. Zr and Ti belong to the same 
column of the Periodic Table, but Zr has a larger effec- 
tive radius and will therefore tend to expand the lattice 
constant of PZT relative to pure PT; from the point of 
view of a Ti ion, this could produce a similar effect as is 
obtained from application of a negative pressure to pure 
PT. 

Here, we explore the strain response as a function of 
both electric field and negative pressure, with and with- 
out the epitaxial constraint. The exercise is carried out 
in three steps. First, we apply a negative pressure in 
the presence of the epitaxial constraint, obtaining the 
strain response as a function of electric field. (In this 
case only the zz component of the stress tensor is rele- 
vant, and we could equally well say that we are applying 
a tensile uniaxial stress along c.) Second, we remove the 
negative pressure and repeat the calculations presented 
in Sec. IIII Al but without epitaxial constraint, so that 
the crystal is free to relax its in-plane lattice constants 
in response to the applied field. Finally, we combine the 
two steps described above and study the behavior as a 
function of electric field and negative hydrostatic pres- 
sure in the absence of any epitaxial constraint. Some 
results of this kind, both at zero and negative pressure, 
have already been presented in Ref. 0. 

Our present results for the application of negative pres- 
sure in the presence of the epitaxial constraint arc shown 
in Fig. [H where the strain plotted on the vertical axis 
is defined as relative to the c lattice constant at zero 
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FIG. 5: Field-induced strain, relative to out-of-plane lat- 
tice constant of zero-field zero-pressure state, for tetrago- 
nal PbTiO-3 under free-stress (0 GPa) and negative-pressure 
( — 1 and —2 GPa) boundary conditions, calculated with 
LAUTREC. Top curves: out-of-plane strain. Bottom curves: 
in-plane strain. 



pressure and zero field. In this and subsequent figures 
we present calculations performed with the LAUTREC 
package. The zero-pressure curve duplicates the data 
presented in Fig. [TJ while the ones at —1 and —2 show 
the enhancement in the strain caused by the negative 
pressure, which shifts the curves by about 2% per GPa. 
Otherwise they look rather similar, except that there is 
slightly more nonlinearity in the curves as the pressure 
becomes more negative. 

The removal of the epitaxial constraint causes a much 
more drastic change in behavior, however, as shown by 
the square symbols (red curves) in Fig. Both in-plane 
and out-of-plane strains are defined relative to the c lat- 
tice constant at zero field. As expected, the enhanced 
polarization and enhanced tetragonality induced by the 
field cause the out-of-plane lattice constant to grow, and 
the in-plane one to shrink, with increasing field. How- 
ever, this variation shows an anomalous behavior: the 
change in strain starts out slowly at lower fields, accel- 
erates and occurs most rapidly at a characteristic field 
£ anom ~ 550 MV/m, and then slows again at higher 
fields. This anomalous behavior is most evident in the 
curve for the out-of-plane strain, but is also visible for the 
in-plane strain. The anomaly also shows up clearly in the 
behaviors of the piezoelectric coefficient and internal dis- 
placements, presented in Fig. EJa) and (b) respectively. 

Finally, the remaining curves in Fig. O show how this 
behavior evolves as a negative pressure is applied. The 
anomaly becomes more pronounced as the pressure be- 
comes more negative, with -E a nom shifting to lower fields 
of about 375 and 150 MV/m at -1 and -2 GPa re- 
spectively. (This behavior was also presented in Ref. HI 
in terms of reduced electrical coordinates.) Very large 
strains (with c/a ratios approaching 1.25) occur at the 
largest fields and strongest negative pressures considered. 
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FIG. 6: Variation of properties of tetragonal PbTiOa with 
electric field under free-stress boundary conditions, calculated 
with LAUTREC. (a) Piezoelectric coefficient ££33. (b) Out-of- 
plane displacements of ions from zero-field positions. 



The behaviors are similar to those already reported in 
Ref. as a function of negative pressure alone, but here 
we find that the application of an electric field works 
cooperatively with the negative pressure to induce the 
crossover into the supertetragonal state. 5 

Returning to the zero-pressure results, we can now 
compare these more directly with previous works in 
which the epitaxy constraint was absent. As was shown 
in Fig.[6ja), we find that ^33 starts at about 70pm/V at 
zero field, increases up to about 140pm/V near -E a nom = 
550 MV/m, and then falls again at higher fields. The 
results of Ref. [5j instead show d 33 starting at 39pm/V 
and decreasing monotonically as the field is increased. 
The discrepancy with respect to our results may be at- 
tributed to the more approximate methods used in Ref. 5, 
where the electric field was applied in an approximation 
in which extra forces appear on the atoms in proportion 
to the Born effective charges computed at zero field 
The Landau-Devonshire results of Ref. [H also show ^33 
decreasing monotonically with field, but starting from a 
zero-field value of 79pm/V that is closer to ours. The 
comparison between their work and ours deserve further 
comment, as will be provided in the next section. 



C. Discussion 

In this section we discuss some of the factors that need 
to be considered when comparing our theoretical calcu- 
lations with previous theory and experiment. We focus 
first on three important considerations, namely, the elas- 



tic boundary conditions (epitaxial vs. free-stress), the ref- 
erence in-plane lattice constant, and the Zr content. We 
then discuss several other factors that may also play a 
significant role. Some comparisons have already been 
made along these lines in Sec. IIII A[ but we concentrate 
here on possible mechanisms and physical explanations 
for the observed trends. 



1. Reference in-plane lattice constant 

When applying the epitaxial constraint, we have a 
choice of which in-plane lattice constant a to use for the 
comparison. As explained earlier, most of our results are 
reported for a value of a that is 1.14% smaller than the 
equilibrium lattice constant of PT, to facilitate compar- 
ison with the experiment of Grigoriev et alA on PZT 
20/80. If instead we repeat the calculations by fixing 
a to be that of tetragonal PT at its theoretical equilib- 
rium lattice constant, as mentioned earlier near the end 
of Sec. IIII At we find a zero-field ^33 of about 50pm/V, 
to be compared with the value of about 68pm/V ob- 
tained at the smaller lattice constant with LAUTREC 
code. This trend is understandable since one may expect 
a larger e?33 in a more tetragonal material. It is also con- 
sistent with the trend shown in Fig. 4(a) of Chen et at, 
where (I33 increases from 39 to 58 pm/V when the epitaxy 
constraint is made 1% more compressive. 



2. Epitaxial vs. free-stress boundary conditions 

Our work reinforces the expectation that the choice of 
epitaxial vs. free-stress elastic boundary conditions plays 
a crucial role in the piezoelectric response. This is true 
already at zero field, but the effect becomes enormous for 
elevated fields. 

At zero field, one expects the free-stress to exceed 
the epitaxial one, other things being equal. This can 
be understood as follows. If we start from the relaxed 
tetragonal ferroelectric PT crystal and apply an enhanc- 
ing electric field along z with the in-plane a fixed, we can 
decompose the response according to a two-step process. 
First, we apply the field while allowing a and c to relax. 
We expect the degree of tetragonality to increase along 
with P, so that a shrinks while c grows, as confirmed 
by Fig. [5l Second, we enlarge a back to its zero-field 
value, again while allowing c to relax; since the Poisson 
ratio is positive, this should cause c to shrink. Thus, 
c?33 should be smaller in the epi-constrained case. This 
is true in the work of Chen et al£>, where for example 
the zero-field d^ decreases from 79 to 31 pm/V for pure 
PT, and from 87 to 45 pm/V for PZT 20/80, when going 
from free-stress to epi-constrained boundary conditions. 
It is also confirmed in our calculations on pure PT; we 
find that the zero-field d^ decreases from 71 to 50 pm/V 
(LAUTREC values) in going from the free-stress to the 
epi-constrained case. 
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As indicated earlier, we find a drastic difference in the 
non-linear response when we remove the epitaxy con- 
straint, with the smooth decrease of ^33 in Fig. EJa) re- 
placed by the rapid increase and peak around 550MV/m 
in Fig. EJa). This peak is clearly a signature of the 
crossover into the supertetragonal state. 



the influence of the perovskite-electrode interface o 31 i 32 
will also play an important role. We expect that these 
issues will be clarified as improved methods of sample 
preparation become available, and as systematic studies 
are done to see how the piezoelectric responses depend 
on film thickness and other properties. 



3. Zr content 

Of course, it is important to keep track of the differ- 
ences between pure PT and PZT with different Zr con- 
centrations, as we have tried to do consistently above. As 
is well known, PZT is often preferred as a more practical 
material for experimental purposes and for applications 
because of reduced leakage currents and other beneficial 
properties. The work of Chen et al£> reports ^33 increas- 
ing from 79 to 87pm/V in going from pure PT to PZT 
20/80 in the free-stress case, and from 31 to 39pm/V in 
the epi-constrained case at zero misfit. If this trend is 
correct, it is in the wrong direction to explain the differ- 
ence between our theory and the experiment of Grigoriev 
et aliZ, since our CZ33 values for pure PT are larger than 
theirs for PZT 20/80. Other factors may be responsible 
for this discrepancy, as discussed in the next subsection. 
Ideally it would be advantageous if the pulsed-field ex- 
periments could be carried out on a series of samples 
of varying Zr content so that an extrapolation could be 
made to the pure-PZ case, allowing a more direct com- 
parison with theory. 

4- Other factors 

Here we comment briefly on a number of other factors 
that might play a role when comparing our results with 
experiment. Our theory is purely a zero-temperature the- 
ory, and also completely neglects the effects of disorder 
arising from Zr configurations in PZT, or from defects 
such as oxygen vacancies even in pure PT. Both dis- 
order and thermal fluctuations may tend to cause the 
local polarization direction to fluctuate about the global 
tetragonal axis, and it is well known that polarization ro- 
tation can generate large electromechanical responses in 
this class of materials—^ At larger length scales, it is also 
possible that the samples might not be in a perfect single- 
domain state, with the incomplete switching being caused 
by defects, roughness, or other imperfections in the thin- 
film samples. In such a case, the intrinsic piezoelectric 
response would be underestimated, while extrinsic contri- 
butions associated with domain-wall motion might also 
be present. For ultrathin films, depolarization effects and 



IV. SUMMARY 

In summary, we have studied the piezoelectric response 
of PbTiOs to a polarization-enhancing electric field ap- 
plied along the tetragonal axis under several kinds of 
mechanical boundary conditions. In the epitaxially con- 
strained regime we find hardly any saturation of the 
piezoelectric coefficient (^33 up to a field of 500MV/m, in 
agreement with recent experimental measurements and 
in contrast with the predictions of Landau-theory expan- 
sions. With the removal of the epitaxial constraint we 
find a remarkable non-linear effect, with ^33 rising to 
twice its zero-field value at ~550MV/m and then de- 
creasing again for higher applied fields. 

The comparison between ABINIT and LAUTREC 
shows fairly good agreement. Some quantitative discrep- 
ancies do exist (e.g., in equilibrium lattice constants); 
these can be traced to the sensitivity of PbTiOs to the 
choice of pseudopotentials. Nevertheless, the trends in 
the two calculations are very similar. Moreover, exten- 
sive tests show that the dissimilarities in the finite-field 
techniques (constrained- E with Berry phase in ABINIT 
vs. constrained-!? with Wannier functions in LAUTREC) 
have little or no impact on the calculated properties, once 
the subtleties in the treatment of the electrical and strain 
variables have been properly accounted for. 

Quantitative comparisons with experiment still present 
a formidable challenge, in part because it remains diffi- 
cult to characterize the precise sample conditions under- 
lying a given electrical measurement. However, signifi- 
cant progress has now been made on the theoretical side, 
and experimental control of film properties continues to 
improve. Thus, we hope that direct comparisons between 
theory and experiment regarding nonlinear piezoelectric 
behavior will provide an increasingly fruitful avenue for 
future investigations. 
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